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I. Introduction 

The advantages, disadvantages, and possible 
applications of H2*02 cambustion driven MHD power 
generators have been discussed previously (refs. 
1-4). Of prime Interest to the NASA Lewis Hz -02 
MHD program Is the development ol llghtwel^t, 
multi -megawatt MHD power supplies. The major dif- 
ficulty in their oevelopaent Is their small size 
which makes surface losses (thermal and electrical) 
a significant fraction of the power generated. To 
overcome these losses, these generators must be 
operated at high Mach numbers and high magnetic 
field. Therefore, our Initial experimental pro- 
gram Is directed toward evaluating the performsmcu 
of H^-Og generators at high voltage, l.e., high 
Mach number (2), and hlgli magnetic field (7 Tesla). 

Two facilities are being used In the NASA 
Lewis Hg-Og experimental progrsun. The primary fa- 
cility is the MID power generation facility which 
was previously described In reference 1. The main 
feature of this facility is a liquid neon-cooled 
cryomagnet capable of a 7 Tesla peak field. At 
present, constjructlon of this facility has been 
completed and systems checkout Is In progress. 

The second facility is a rocK’" test facility 
for testing MHD hardware and performing diagnostic 
tests on the alkali-metal seeded combustion plasma 
and MHD duct. This facility has no magnet end, 
therefore, diagnostic measurements are made with 
applied electric fields. This facility is shown in 
operation in figure 1. Shown are the injector- 
combustion chamber -nozzle and a prototype MHD Hall 
generator. At present, over 1000 test runs for 
times ranging from a few seconds to 20 seconds have 
been performed in this facility. These tests were 
performed to develop Intcrelectrode Insulating ma- 
terials, to measure heat losses, and to measure the 
fluid dynamic and electrical properties of MHD 
ducts ana the cesium seeded Hg-Og working fluid. 

II. Combustion Chamber-Nozzle Heat Loss Teste 

An initial concern with the water-cooled cop- 
per combustion chambers (including supersonic noz- 
zle) was that these chambers were designed to op- 
erated steady state with wall temperatures some- 
what beyond the boiling point of water, but con- 
siderably below the conde.nsat ion point for CSOH. 
These chambers therefore have two undesirable fea- 
tures due to the low wall temperature; possible 
high heat loss and condensation of CSOH which could 
then short out the MHD channel by streaming down 
the channel as a liquid. To alleviate these prob- 
lems, a second set of combustion chamber -nozzles 
was constructed with a 0. 006-inch graded mixture of 
zirconia 'niohrome flame sprayed on their inner sur- 
face. This provides an allowable surface tempera- 
ture well above the condensation temperature of 
CSOH. 


singly the heat transfer to the hot wall ceramic- 
coated chambers is higher than to the cold wall cop- 
per chambers. Optical pyroaeter measurements in the 
supersonic portion of the combustion chamber-nozzle 
verified that the ceramic coating did produce a sur- 
face temperature of approximately 1500” K in this 
region. 

To gain some understanding of these results, a 
simplified heat transfer calculation based upon un- 
accelerated flows was carried out. This simplifica- 
tion leads to an overestimate of the heat flux in 
the acceleratlng-deaccelerating region around the 
throat of the nozzle. However, since the surface 
area of this region is relatively small compared to 
that of the total surface area of the combustion 
cbo.r^r-nozzle, the effect upon the net heat trans- 
fer should be small. (For a discussion of this 
point see ref. 5.) If one assumes fully turbulent 
flow throughout the combustion chamber-nozzle, the 
heat flux, q, is given by 

q - 0.026 I (Ke)°-®(Pr)°-^(T^^ - T^) (1) 

where: k is the thermal conductivity 

D is the duct diameter at the point at 
which the heat flux is evaluated 
Re is the Reynolds .number 
Pr is the Prandtl .number 
Tad is the adiabatic temperature 

3 T + (Pr)l/3(Ts - T) (2) 

T is the static temperature of the 
flowing gas 

T3 is the stag.nation temperature of the 
flowing gas 

is the ccmoustor-nozzle wall tempera- 
ture 

and k, Pe, and Pr are evaluated at the reference 
temperature condition, 

Tf = T + 0.5(T„ - T) + 0.22(T^d - T) (3) 

The above equations are those of reference 5 
for vjiaccele rated flows except they are for a heat 
tremsfer coefficient based upon temperature rather 
than enthalpy. This was done so that tabular data 
from reference 6 could be used to evaluate viscos- 
ity, thermal conductivity, and heat capacity. 

These coefficients were evaluated assuming chemical 
equilibrium through ^'he boundary laiyer. 

The heat flux from the gas (eq. (1)) is then 
equal to the heat flux throxigh the chajiiber wall 
given by 



The net heat loss from these two designs was 
determined by measuring the cooling water flow rate 
and temperature difference from inlet to outlet. 

The measured values are shown in figure 2 as a 
function of combustion chamber pressure. Surpri- 




where: Tj^ is the temperature of the cooling water 

and '^re the thickness/ 

1 V /coat ' 'cu 

I thermal .-onductivity for the ceramic 

coating and copper wall, respectively 
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Is the filB coefficient for the cooling 
water 

The properties of the flowing gas are dater- 
Blned by Integrating the pressure and temperature 
equations for the flowing gaa through the contour 
of the ccnbustlon chaoher-nossle. These equations 
Include heat transfer (eq. ().)), friction (smooth 
pipe) and equilibrium chemistry (ref. 7). The 
problem of the singular point at Mach > 1 la re- 
moved by integrating subsonlcally until (1 - Mach) 

< 0.001 and then Jumping across to Mach al.OOl. 

The Integration technique Is Runge-nitta-Merson 
(ref. e) whose variable step size allows the treat- 
loent of the singularity to be prograaned Into the 
logic. 

The heat transfer at each station of the Inte- 
gration Is determined by iteratively solving equa- 
tions (1) to (4) for the unknown temperatures 
Tv, Tj., and the heat flux q. Note that X an,. 

T, sire given by the gas dynamic calculation and 
cooling water temperature is taken to be room 

teaqwratiu'e . 

The results of this calculation are shown by 
the solid lines In figure 2. The agreement with 
experiment is good and Indeed the theory predicts 
that the total heat transfer to the hot wall 
ceramic coated combustors Is greater than that tc 
the cold wall copper combustors. Examination of 
the results of the theoretical calculations indi- 
cated that the reason for this apparent anomaly is 
that while the factor Tj^^ - Tv (eq. (1)) decreases 
with Increasing wall temperature, i.e., for the 
ceramic coating, this decrease is more than off- 
set by the increase in the thermal conductivity 
resulting from the increase in Tp (eq. (3)). 

III. fluid I!y n«I-^ r TAB-hs 

Another concern was with the fluid dynamic 
performance of the nozzle. The design was adapted 
from a previous rocket program and Involved a com- 
promise between achieving a low weight and per- 
fectly parallel flow. Careful pressure probing at 

0.635 cm (0.25 in.) intervale from the end of the 
nozzle to a point approximately 8 cm downstream In- 
dicated the presence of a pressure disturbance ap- 
proximately 5 cm downstream of the nozzle exit. 

The data a^ shown in figure 3. The symbols un the 
figure correspond to the following channel configu- 
rations. 

1. C 12-0.635 cm (0.25 In.) electrodes with 
pressure taps on electrodes 1, 3, 5, 7, 9, and 11 
as numbered from the nozzle 

2. X 2.54 cm (1 In.) spacer between nozzle 
and the above channel 

3. O channel as In 1 except turned end for 

end 

4. + 30-0.127 cn (0.5 In.) electrodes with 
2.54 cm spacer between nozzle and channel with 
pressure taps on all but electrodes 10, 11, and 20 

To understand the nature of this disturbance, 
the experimental conditions were ansulyzed using a 
three-dimensional supersonic nozzle analysis 
(ref. 9). This Is a modification of the scheme of 
reference 10 which Is an exact solution to the 
chforacterlstic equations and is accurate to second 
order In the grid spacing. While the analysis 


smooths out the sharp changes across a shock wave 
and ignores all dissipative effects, reference 9 
shows that it predicts the location and siagnltude of 
strong compressions fairly well. The theoretical 
prediction (shown in fig. 3) Is seen to agree well 
with the experimental data as to the location and 
magnitude of the first oblique shock. Some recent 
data taker In a MHD duct built for initial tests in 
our mag:. St facility are also shown by the crossed 
points in figure 3. While the data are of insuffi- 
cient resolution to establish the exact structure 
of the downstream pressure disturbances, there is 
evidence of at least quantitative agreement with 
theory. 

Also shown on figure 3 Is a one-dismnslonal 
MHD channel calculation Including friction, heat 
transfer, and chemistry. This calculation is simi- 
lar to that described In the prevloua section with 
the exceptloti that in keeping with the results ob- 
served by others (ref. 11 and references included 
therein), the skin friction coefficient was in- 
creased by a factor of 3 over the smooth wall value 
to compensate for the ro\igh wall struct\ire of the 
segmented MHD duct. Again, the agreement is good. 
Note that the agreement is not with an average 
through cbe data since the distance over which the 
shocks occur Is small coiqpared to the distance be- 
tween shocks. Therefore, the theory should be more 
heavily weighted by the points in the valley rather 
than the peaks of the pressure disturbance (ccBq>ere 
1-D isentroplc auid 3-D nondissipative) . 

While the observed pressure disturbances do not 
appear to cause additional pressing drop over that 
resulting from rou^ walls, their effect iqion the 
performance of the duct in the presence of a magne- 
tic field Is open to question. The pressure dls- 
t’urbances can lead to nonuniformities In the elec- 
trical conductivity and/or Kali parameter which in 
turn may lead to a deterioration In performance of 
the generator In the presence of a magnetic field. 
Such deterioration has been predicted theoretically 
for plane, layered nonuniformities in these param- 
eters (see ref. 11 and references cited therein). 

IV. Electrical Conductivity Tests 

Electrical conductivity tests were performed 
In constant area, heat sink, cylindrical, Hall type 
MHD channels. The electrical conductivity was icez- 
sured In the classical manner (see ref. 12) by ap- 
plying a voltage across the channel from one end 
electrode to the other, measuring the current, and 
usi.Tg the inner electrodes as probes to monitor the 
voltage distribution along the channel. Typical 
voltage profiles are shown In figure 4 as a function 
of time. The channel used in this test consisted of 
12-1.27 cm (0.5 In.) wide electrodes separated by 
synthetic mica Insulators. It is noted from this 
flgu -e that there are voltage drops at both ends of 
the :hannel with a linear region In the central por- 
tion, The voltage drops ajre associated with the 
current flow across the cold boundary layer and then 
into the larger area associated with uniform current 
flow through the channel. Therefore the linear re- 
gion Is representative of the electrical properties 
of the gas and hence was used to Infer the avereige 
electrical conductivity of the gas from 

I 

° ’ A(AV/iX) 




vhert: o ia the electrical conductivity 

K la the croaa sectional area of t> a 
channel 

eH/eX. la the voltage gradient dete.^dned 
from the voltam curvea 
I la the neaaured current 

For the teat data ahown In figure 4, the 
working fluid was a atolchlooetrlc mixture of 
Hg-Og at a ccohuatlon pleasure of 1.034x10^ R/M^ 
(150 psla). The Mach 2 flow was seeded with a 75^ 
solution of CSQH dissolved in water. The solution 
was sprayed Into the oxygen feed line. For this 
test the following data were obtained. 


Time (sec) 

o(MH0/M) 

I (Amp) 

AV/AX (volts/M) 

2.08 

12.62 

1.762 

72.4 

3.08 

13.16 

2.037 

80.3 

4.08 

12.89 

2.112 

85.0 


where aV/AX Is the vedue in tl e linear region 
(electrodes 3-8). 


It is noted that the value of the gas conduc- 
tivity does not change appreciably during the run. 
However, the measured current varies from 1.762 to 
2.112 amps due to the Increase In AV/aX. This Is 
a result of the decrease In the electrode voltage 
drop at the ends of the channel. This decrease In 
voltage drop is the result of the Increasing wall 
temperature of the heat sink channel with run time 
and the consequent increased temperature of the 
boundary layer. 

In figure 5 the results of several conduc- 
tivity measurements are given. The solid curves 
are theoretical values calculated usl.cg the chem- 
istry of reference 7 to calculate the free elec- 
tron density taking into account a heat loss in 
the combustor-nozzle and a 1^ heat loss In the 
channel. The calculated values are the local val- 
ues at the midpoint of the channel. The only im- 
portemt contributions to the collision frequency 
ai-e from HgO and cesium compo’i.'ids. The HgO cross 
section was taken from reference 13. Cesium and 
its compounds w«>ra given the cross section from 
reference l4. Cross section values were assigned 
to n-i 1 other species, but due to their unlsq>or- 
tance and number (31) they will not be referenced. 

The seed fractions over which conductivity 
data has so far been taken Is limited by the size 
of the seed flow meter. Therefore as noted In 
figure 5, the range changes as we Increase the 
combustion pressure, l.e., H2-0g mass flow. How- 
ever, the following conclusions can be drawn. For 
a combust lor pressure of 2.069x10® N/^^ (300 psla) 
the experimental data are In good agreement with 
theory. At a combustion pressure of 1.034x10° 

N/M^ (150 psla) the experimental data agree with 
theory up to a seed fraction of approximately 0.05 
at which point It decreases well before the expec- 
ted maximum at a seed fraction of 0.15. The ex- 
perimental points at a combustion pressure of 
5,17x10® (75 psla) are nearly a factor of 5 below 
thqory and show no dependence on seed fraction. 

The lower than predicted electrical conduc- 
tivity at higher seed fractions and/or lower pres- 
sures is thought to be due to poor atomization of 
the seed as the ratio of seed mass flow/oxygen 
mass flow increases. This effect is well known 


from studies of two-fluid atCBlzatlon, e.g, , ref- 
erence 14, where it Is shown that the droplet size 
increases rapidly as the ratio of liquid to gas 
mass flow Increases. In our experiment these 
Isurger droplets may be of sufficient size that they 
are not vaporized during their residence time In I 
the coBibustor. This effect would be Increasingly 
noticeable at lower pressures as is observed In our 
experiments, 

V. Coneludlng Remarks 

This study has shown that flame sprayed cera- 
mic coatings on combustion chamber and nozzle walls 
permits operation with hot walls, thereby elimi- 
nating possible problems associated with seed con- 
densation. However, under the conditions of our 
experiment the net heat transfer was unexpectedly 
Increased over that of cold copper walls. Further 
work will be concerned with compari.ng the electri- 
cal conductivity of the gas at the exit of these 
nozzles. 

Pressure disturbances associated with oblique 
pressure waves have been measured In our MHS chan- 
nel. These disturbances are due to a rocket-type 
nozzle which was designed to achieve low .eight 
rather than perfectly parallel flow at the exit. A 
new nozzle with parallel flow at the exit has ) een 
designed and will be used In our mag.net fa;ll''cy 
along with the previous nozzles to study t. - effect 
of t' e pressure disturbances upon the performance 
of the MHD generator. 

Electrical conductivity measurements were 
found to be in good eigreement with theory except at 
low ccmbustlon pressures and/or high ratios of 
seed/ oxygen mass flows (seed was injected into oj^- 
gen flow line). The discrepancy is believed to be 
the result of poor atomization of the seed at high 
ratios of seed/ oxygen mass flows which results In 
droplet sizes which do not ccmpletely vaporize dur- 
ing thel'* residence time in the combustor. This 
problem is being analyzed spectrogr. phically. 
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Figure 1. - Rocket Lab test facility. 
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Figure 2. - Combustor- nozzle heat loss versus combus- 
tion pressure. 
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Figure 3. - Static pressure distribution along supersonic nozzle 
section and channel. 






